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Direct observation of optically induced transient structures in graphite using ultrafast 

electron crystallography 

Ramani K. Raman, Yoshie Murooka, Chong-Yu Ruanjj Teng Yang, Savas Berber, and David TomaneltiJ 

Physics and Astronomy Department, Michigan State University, East Lansing, Michigan 48824-2320 

We use ultrafast electron crystallography to study structural changes induced in graphite by a 
femtosecond laser pulse. At moderate fluences of <21 mj/cm^, lattice vibrations are observed 
to thermalize on a time scale of Ri8 ps. At higher fluences approaching the damage threshold, 
lattice vibration amplitudes saturate. Following a marked initial contraction, graphite is driven 
nonthermally into a transient state with sp'^-like character, forming interlayer bonds. Using ab 
initio density functional calculations, we trace the governing mechanism back to electronic structure 
changes following the photo-excitation. 

PACS numbers: 82.53.Mj, 81.05. Uw, 61.14.-x, 65.40.De 
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There is growing interest in displacing atoms in mate- 
rials by photo-excitations [l] . Observations of transient 
structures thus formed offer a glimpse into the transfor- 
mation pathways between different structures. Carbon, 
with its propensity to form a wide range of bonding net- 
works {sp, sp^, sp^) is ideal to study the dynamics of 
bond formation and rupture. Of particular interest is 
the conversion of graphite to diamond [2, 311, which is 
believed to involve the rhombohedral phase of graphite 
as intermediate state [J, |5|. Whereas pioneering ultra- 
fast optical studies of graphite have provided evidence 
for photo-induced melting [y, |7| and generation of coher- 
ent phonons [8|, |9[ by observing changes in the electronic 
properties, direct determination of lattice structural dy- 
namics from optical data proves difficult, especially in the 
far-from-equilibrium regime. X-ray diffraction has been 
successful in observing nonthermal structural changes in 
semiconductors [10[, relying primarily on the integrated 
Bragg intensities. However, direct observation of atomic 
motion in nanostructures with low atomic number, such 
as carbon, has not yet been achieved. Electron diffrac- 
tion, with its five orders of magnitude enhanced scatter- 
ing cross-section and advanced schemes achieving fem- 
tosecond temporal resolution [ll|, [l^, ll3|, llil, offers a 



new window into the realm of photo-excited structural 
dynamics, with sensitivity down to ^1 nm llLll4Lll5l|. 



Here, we report the first direct determination of struc- 
tural changes induced in graphite by a femtosecond laser 
pulse. At moderate ffucnces of <21 mJ/cm^, after a 
short thermalization period of ~8 ps, we can attribute 
a temperature value to the graphitic layers and find the 
interlayer vibration amplitudes to match those reported 
in X-ray and neutron scattering studies [16|. At higher 
fluence values approaching the damage threshold, we ob- 
serve lattice vibration amplitudes to saturate. Following 
a marked initial contraction of the interlayer spacing by 
<6%, graphite is driven nonthermally into a transient 
state with sp'^-like character, forming 1.9 A long inter- 
layer bonds. Using ab initio density functional theory 
(DFT) calculations, we trace the structural changes back 




FIG. 1: Uhrafast Electron Crystallography (UEC) of highly 
oriented pyrolytic graphite, (a) The UEC pump-probe setup, 
(b) The layered structure of graphite, (c) Ground state 
diffraction pattern of graphite, (d) Ground state layer density 
distribution function (LDF), obtained via a Fourier transform 
of the central streak pattern in (c). 



to a nonthermal heating of the electron gas, followed by 
a photo-induced charge separation causing a compressive 
Coulomb stress. 

The experimental setup, shown in Fig. [ija), has been 
described in detail elsewhere [14,]. Freshly cleaved sam- 
ples of highly oriented pyrolytic graphite (HOPG) were 
placed inside an ultrahigh vacuum chamber through a 
load lock system at room temperature. A mode-locked 
Ti-Sapphire laser (p— polarized, 45 fs, 800 nm, 1 kHz) 
generated the pump pulse, which was marginally focused 
(w600 /im diameter) onto the sample at «45° relative 
to the HOPG c-axis. The range of the laser fluence 
0.5 < F < 90 mJ/cm^ was below the in situ determined 
optical damage threshold of «120 mJ/cm^. 

A photo-generated 30 keV electron beam (wavelength 
Ae = 0.069 A), with diameter demagnified to <5 /zm, was 
directed onto the sample at grazing incidence of 2.4° to 
serve as the probe. A translational optical delay stage 
was used to vary the relative delay between the arrival 
of pump and probe pulses at the sample, to observe the 
changes induced by laser irradiation. The electron pulse 
sampled the top 3 — 5 layers of the graphite structure, 



depicted in Fig. [Ub) , and produced oriented molecular 
diffraction patterns in the central streak region, shown in 
Fig. [He). The layer density distribution function (LDF), 
shown in Fig. [U^d) , was obtained from the Fourier anal- 
ysis of the interference pattern. The primary peak at 
3.35 A and less pronounced peaks at 6.7 A and «10 A 
were found to be in good agreement with the layered 
structure of bulk graphite [l6| . The decay of higher or- 
der LDF peaks suggest a probing depth of «1 nm. 

To study the photo-induced structure dynamics, we 
first examine the near-equilibrium regime at low fluences 
of 0.5 mJ/cm^<F<21 mJ/cm^. Random atomic dis- 
placements diminish the intensity of the diffraction spec- 
tra by a Debye- Waller factor e~^*^, where M = —s^u^ /A, 
V? is the mean-square atomic displacement perpendicular 
to the reflecting planes, and s = (47r/Ae) sin(0/2) is the 
momentum transfer associated with the maxima located 
at the scattering angle 9. Thus, the change in the mean- 
square atomic displacement Am^, measured relative to 
the unperturbed state at negative times, determines the 
rise in lattice temperature and can be calculated from the 
diffraction intensities as ln(/s(i)//s(i < 0)) — — s^Au^/4. 
In Ultrafast Electron Crystallography (UEC), we moni- 
tor the integrated intensity of the (0,0,6), (0,0,8) and the 
(0,0,10) maxima arising from the interference between 
the graphitic layers, shown in Fig.[2Ka). The intensity of 
all 3 maxima drops within a timescale of 8±1 ps, indica- 
tive of increased thermal motion of the graphitic planes. 
As we probe reflections from the basal planes, the ob- 
served intensity drops mainly due to the increased out- 
of-plane displacement of the atoms. Recent optical stud- 
ies [Sl |9| have associated near-infrared optical excitations 
in graphite with a femtosecond generation of coherent 
phonons with E2g symmetry (inter layer shearing mode). 
Hence, this 8 ps timescale of inter layer thermal excitation 
is a direct measure of the phonon-phonon interactions in 
HOPG, and is in good agreement with reports of hot 
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phonon relaxation times in graphite 

Quantitative analysis of the intensity drops in Fig. [2ta) 
indicates that the mean square atomic displacement lS.v? 
increases linearly with the applied fluence in the near- 
equilibrium regime, as seen in Fig. [2Kb). Comparing our 
data to a theoretical model of temperature dependent 
Au^, benchmarked with X-Ray and neutron scattering 
data [16| , we estimate a temperature rise of 950 K at F = 
21 mJ/cm^ in Fig. [^b). This is also in good agreement 
with temperature extracted from optical studies using 
heat capacity and absorbance [71]. 

As we enter the far-from-equilibrium regime by increas- 
ing the excitation fluence beyond 21 mJ/cm^, the am- 
plitude of atomic vibrations saturates and no longer in- 
creases with increasing fluence, as seen in Fig.[2Kb). This 
is rather surprising, for one would expect a monotonic in- 
crease in the thermal motion as graphite approaches its 
melting threshold at 130 mJ/cm^ 7j. Furthermore, the 
(0,0,6) maximum is found to shift by almost 0.34 A~^ 
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FIG. 2: (a) Decay of the diffraction intensity, (b) Change 
in mean-squared atomic displacements perpendicular to the 
graphite planes, estimated from the Debye- Waller analysis of 
the intensity drops in (a), (c) Effective surface potential Vs at 
the HOPG surface as a function of the applied fluence F. (d) 
Decay of the effective surface potential with fit to the drift- 
diffusion recombination model (see text), (e) Contraction of 
the interlayer spacing Ad (left axis) and the corresponding 
increase in the effective surface potential (right axis) in the 
far-from equilibrium regime. 



at i^ = 40 mJ/cm^, which would correspond to an inter- 
planar expansion of 6%. Such a large peak shift cannot 
arise from structural changes alone. 

We believe that the key effect of the optical excita- 
tion is the repopulation of electronic states, causin g re - 
bonding and redistribution of charge at the surface [l8j. 
This charge redistribution gives rise to a Coulomb field 
near the graphite surface, which induces a collective shift 
of the electron diffraction pattern. We measure this 
Coulomb field here directly by modeling the 'refraction' 
effect commonly associated with the existence of a sur- 
face potential at the vacuum-material interface [18|. We 
deduce the effective surface potential Vs as a function of 
the applied fluence and find it to increase linearly up to 
FwSO mJ/cm^, as shown in Fig. [D^c). This linear trend 
indicates that the surface charging is mainly driven by 
the non-equilibrium diffusion of photogenerated carriers 
rather than multiphoton ionization. The decay of this 
transient surface potential resulting from the recombina- 
tion of space charge follows a power-law decay with the 
exponent —1, as shown in Fig. 2(d) and predicted by the 
drift-diffusion recombination model [18| . 

To examine the effect of optical excitations on the 
atomic structure at these elevated fluences, we monitored 
the interlayer separation following the photo-excitation 
by a time-resolved LDF analysis, after having accounted 
for the surface potential effects. In contrast to the ther- 
mal expansion at low fluences, we observed the interlayer 
separation to contract rapidly, as seen in Fig. [^Je) . The 
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FIG. 3: Signature of sp^-like bonding in HOPG. (a) 
Molecular interference pattern M{s) and the corresponding 
LDF curves at selected time stances following strong photo- 
excitation. The transient peak at 1.9 A in the LDF indicates 
formation of an interlayer bond, (b) Ratio of peak intensities 
obtained from the diffraction patterns. The enhancement of 
the (0, 0, 12) with respect to the (0, 0, 10) peak cannot be ex- 
plained by thermal vibrations and indicates a transient struc- 
ture consistent with sp bonding, (c) Ratio of peak intensities 
at 7? = 1.9 A and 3.35 A in the LDF. 



change in the interlayer spacing Ad is found to correlate 
well with the rise and fall of Vg a.t F — 77 mJ/cm^, with 
a maximum value T4 ~ 12 V corresponding to an inter- 
layer contraction of « 6%. The potential rise Vs sampled 
within the electron probe depth of 1 nm yields an inter- 
nal field of -E ~ 1.2 V/ A, which causes Coulomb stress. 
Using this value of E and Cr = 10 for graphite, we esti- 
mate a maximum energy density of Ue = ereoE'^/2 = 
0.2 eV/atom, close to the wO.3 eV/atom value based 
on DFT for the activation barrier in the graphite-to- 
diamond transformation [3|, l4| . 

To further elucidate the structural evolution, we exam- 
ined closely the diffraction spectra. Figure[3ja) shows the 
time evolution of the molecular interference curves M(s) 
and the corresponding LDF curves at F = 77 mJ/cm^, 
just below the damage threshold. Starting at ~6 ps, 
as Vs rises beyond 4 V, we observe the onset of peak 
shouldering and broadening in the M{s) curve and ap- 
pearance of a new peak at i?Ril.9 A in the LDF curve. 
Concurrent with this, we also observe a drastic change 
of relative intensities between the even and odd peaks in 
the M{s) curves. In graphite with the interlayer spacing 
of do = Co/2 = 3.35 A, the even Bragg peaks n = 4, 6 are 
labeled (0, 0, 8) and (0, 0, 12), and the odd peaks n = 3, 5 
are labeled (0,0,6) and (0,0,10), respectively. The ob- 
served higher intensity of the (0, 0, 12) peak with respect 
to the (0, 0, 10) peak in Fig. [3jb) cannot be explained 
by thermal vibrations and suggests a significant struc- 
tural change. Consistent with the fact that only the even 
diffraction peaks should be present in the structure of di- 
amond, and also with the observation of a well-defined 
interlayer bond at i?Ril.9 A, we conclude that the tran- 



sient structure is consistent with sp^ bonding. We find 
the transient sp'^-like structure, associated with the rel- 
ative intensity of the i?»1.9 A peak in the LDF curve, 
to reach a maximum fraction of 47% at 14 ps, as seen in 
Fig. [3] (c). Similar sp^ — sp^ hybrid structures have also 
been identified in pressure-driven pathways of graphite- 
to-diamond transition[19|. By 45 ps, the structure re- 
covers its sp^ character, but remains hot even after 3 ns 
with an average value Au^ = 0.033 A^, corresponding to 
a lattice temperature T«1000 K. 

To clarify the origin of the structural changes, we stud- 
ied the effect of laser pulses on the electronic struc- 
ture and bonding in hexagonal graphite using ab ini- 
tio DFT calculations in the local density approxima- 
tion (LDA). While much of the basic physics underlying 
photo-induced structural changes can be understood in 
the bulk structure, we used graphite slabs representing 
the surface for quantitative predictions. The total energy 
was determined using the ABINIT plane-wave code [20| 
with a 64 Ry energy cutoff, TrouUicr-Martins pseudopo- 
tentials, and the Ceperley- Alder form of the exchange- 
correlation functional. The Brillouin zone of the 4— atom 
bulk unit cell was sampled using a fine mesh of 24 x 24 x 12 
fc-points that included the K — H line at the Brillouin 
zone edge, close to the Fermi surface. This approach cor- 
rectly reproduced the observed |16j in-layer bond length 
dec = 1-42 A and the interlayer spacing do — 3.34 A 
in the bulk system at T = 0, shown schematically in 
Fig. m^a). The electronic density of states of graphite 
near the Fermi level is shown in Fig.[4l(b), and the total 
charge density is depicted in Fig. SJc) in a plane normal 
to the layers. In the following, we will study the effect of 
two types of photo-induced electronic excitations on the 
structure. 

First, we will consider the thermalization of the ini- 
tial nonequilibrium population of electronic states in the 
laser-irradiated target over a sub-picosecond time scale 
21 1 to that of a very hot electron gas [22|. The elec- 
tronic temperature (Tg) is probably lower than the lim- 
iting value kBTe^hv = 1.55 eV and depends on the laser 
fluence. For the sake of illustration, we compare the 
Fermi-Dirac distribution in graphite at ksT^ — and 
ksTe = 1.0 eV in Fig. [41(b), indicating the electronic 
excitations in the tt — tt* manifold of the hot electron 
gas. Populating initially empty conduction states by va- 
lence electrons leads to a change in the total electron 
density Ap{r) — p(r; fc^Te) — p{r;0), which is depicted 
in Fig.m^d) for fc^Te = 1.0 eV. 

As stipulated by DFT, changes in the charge density 
modify the force field in the system. Our results in 
Fig-IHd) suggest an increased population of C2pz orbitals 
that may hybridize to ppa bonding states connecting 
neighboring layers, increasing their attraction. The de- 
population of in-layer bonding states, on the other hand, 
should cause an in-layer expansion. To obtain a quanti- 
tative estimate of photo-induced structural changes, we 
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FIG. 4: (Color online) Electronic structure changes during 
photo-excitations in graphite, (a) Schematic view of the struc- 
ture, (b) Electronic density of states (solid line) and the 
Fermi-Dirac distribution of graphite at T = (dashed line) 
and ksT = 1.0 eV (dotted line), (c) Total pseudo-charge den- 
sity p(r) and (d) change in the pseudo-charge density Ap(r) 
due to an electronic temperature increase to fcsT = 1.0 eV. 
The plane used in (c) and (d) intersects the graphene layers 
along the dashed lines in (a). 



performed a set of global structure optimization calcula- 
tions of graphite with electrons subject to effective tem- 
peratures in the range eV< ksT < 1.55 eV. We used a 
stringent convergence criterion, requiring that all stress 
components lie below 5.0x10"^ Ha/ajj and that no force 
exceeds 5.0x10^''' Ha/ap. In bulk graphite, we indeed 
found that nonzero electronic temperature leads to a 
maximum interlayer contraction Ad/ do — —1% and in- 
layer expansion Arcc/fcc — +1% for fcsTeRil.O eV. In 
the slab geometry, we obtained a slightly larger interlayer 
contraction of up to Ad/da = —1.5%. 

To further take into account the effect of Coulomb 
stress induced by the laser pulse, we represent the charge 
separation in the electronic ground state, corresponding 
to the observed internal field of Ril.2 V/A. To model 
this system, we immersed three-layer graphite slabs, sep- 
arated by 30 A, in a uniform electric field and adjusted 
its strength to reproduced the observed field value. With 
the polarized charge distribution frozen in, the external 
electric field was switched off. Then, the internal field 
was only due to the polarized charge distribution. The 
slab calculations were performed using the SIESTA code 
[23J with double-polarized triple-^ local basis and a fine 
24x24 fc-point mesh. The charge density was obtained 
on a real-space grid with a mesh cutoff energy of 250 Ry. 
From total energy calculations at different interlayer sep- 
arations, we found that charge density redistribution as- 
sociated with the internal field may further reduce the 
interlayer separation by 2 — 3%. Combined with the con- 



traction induced by the initial non-equilibrium electron 
heating, we thus can explain the observed contraction of 
the topmost interlayer separation by up to «5%. 

In conclusion, a nonthermal pathway of photo- 
induced structural changes in graphite has been ob- 
served through multidimensional crystallographic deter- 
mination in space, time and energy. Beyond a threshold 
fiuence, a transient sp^-\ike structure with well-defined 
interlayer bonds emerges. The modified force field in the 
excited state and the Coulomb stress are the main forces 
driving this structural change. 
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